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We study the differential conductance in the Kondo regime of a quantum dot coupled to multiple 
leads. When the bias is applied symmetrically on two of the leads (V and —V, as usual in experi- 
ments), while the others are grounded, the conductance through the biased leads always shows the 
expected enhancement at zero bias. However, under asymmetrically applied bias (V and AV, with 
A > ) , a suppression - dip - appears in the differential conductance if the asymmetry coefficient A is 
beyond a given threshold Ao = \fl + r determined by the ratio r of the dot-leads couplings. This is a 
recipe to determine experimentally this ratio which is important for the quantum-dot devices. This 
finding is a direct result of the Keldysh transport formalism. For the illustration we use a many-lead 
Anderson Hamiltonian , the Green functions being calculated in the Lacroix approximation, which 
is generalized to the case of nonequilibrium. 

I. INTRODUCTION 

The observation of the Kondo effect in quantum dots allows a direct insight on this interesting many-body effect, 
whose experimental fingerprint is the enhanced differential conductance at zero bias which is a specific type of zero 
bias anomaly. The great advantage over the bulk Kondo effect is of course the possibility to control the parameters 
like for instance the coupling with the leads. Next, it was shown that even more information can be extracted if the 
dot is connected to three leads. The three lead-Kondo problem was theoretically considered before in, e.g., [H 0, IH 
and experimentally realized by Leturcq et. al @, 0). The third lead was used to read the non-equilibrium density of 
states (DOS) induced by the other two leads. Shah and Rosch Q analyzed the influence of the coupling and bias on 
the decoherence rate. 

In this paper, we also study the many-lead quantum dots in Kondo regime, but the focus is different : we are 
interested to find out the conditions under which the differential conductance can show a suppression at zero bias 
instead of the usual enhancement. There are well-known examples of transport properties influenced by dips in the 
spectral function produced by Coulomb interactions. The first one was the Coulomb pseudo-gap in the impurity 
conduction of doped semiconductors resulting in the T~ 1 / 2 dependence of the hopping conductivity Then, a zero 
bias anomaly manifested as a dip in tunneling DOS was studied in tunnel junctions However, what is conspicuous 
for the three-lead quantum dots (see sketch in Fig. 2), is that a peak-dip crossover may appear in the differential 
conductance, while keeping the dot permanently in Kondo regime. As we show in the next section, this occurs 
when the bias is applied asymmetrically on the left and right leads, and the crossover is triggered by the degree of 
asymmetry. For the sake of definiteness, let us assume that for some voltage configuration (see sketch in Fig. 2), the 
left lead receives electrons from the central lead and pumps electrons into the right lead, the difference of the two 
contributions being the total current. For the conductances, one can write Gl = Glc + Glr = Glc ~ \Glr\- Both 
Glc an d \Glr\ must decay by applying a voltage V (as the voltage destroys the Kondo resonance and reduces the 
electronic tunneling between leads), but not necessarily with the same derivative ; then, if dGhc/dV > d\GLR.\/dV 
the conductance Gl will increase with V , giving rise to a dip about V = 0. 

It will also be shown that the peak-dip crossover can be used to determine the ratio between the couplings of 
different leads to the dot, and therefore it is of potential experimental use. 

The main conclusions of the paper are extracted from the current formula given by the Keldysh formalism (Eq.l) 
by taking advantage of a simple connection between the convexity of the conductance dl/dV and the convexity of 
the density of states. 

Recently, in the first two-channel Kondo experiment Q , a dip generated by the competition between channels was 
observed in the differential conductance, even in the case of symmetrically applied bias. One may note that in our 
calculations, as in the mentioned (theoretical and experimental) papers referring to the 3-lead problem 0,0,0,0); 
the set-up and physics are different, and therefore the results are outside the multi-channel Kondo phenomenology. 

For the illustration of the results, we shall give also some numerical calculations using the equation of motion 
solution for the Anderson Hamiltonian Q , generalized to the multi-lead and non-equilibrium case. 
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II. PEAK-DIP CROSSOVER IN DIFFERENTIAL CONDUCTANCE: ANALYTICAL ANALYSIS 

In order to calculate the electron current through the lead a one can use the Keldysh formalism, and obtain Q : 
I a = / ImG dd {u) ^2 A P _ Ma) - /(^ - Hp)] du, ^ I a = 0, (1) 

where A Q is the coupling of the lead a to the dot, A = A Q , and f(uj — /i Q ) is the Fermi distribution in the 
a lead. G dd is the retarded Green function of the dot in the presence of the leads. Wc remind that, usually, the 
Keldysh formalism needs also the "lesser" Green function , but when the leads are all connected in the same site, the 
lesser function can be eliminated (except for the case of magnetic leads [HI). It was recently proven that the same 
simplification occurs for a symmetric many-site system [ill ] . 

In this paper we are interested in the behavior of the conductance close to equilibrium, for the case when two of the 
leads are biased and the others are grounded. In this case, all the grounded leads have the same chemical potential 
(set to zero) and act like one single reservoir; this reservoir will be called "central lead" (with the chemical potential 
Vc = 0). If we apply the potential Vl — V on the left lead " L" and the potential V R = XV on the right lead "i?", 
then the differential conductance through the lead "L" can be calculated from Eq.l : 

G L (V) = ^ = ^11 [(A* + A c )I mG ^) d ^^ - XA R I mGdd n dfiUJ duj XeV) ]d.. (2) 

Notice that only the derivative of the Fermi function was considered while the derivative of ImG dd versus bias was 
neglected, which was proven to be a very good approximation in Q. At low temperatures, when the derivatives of 
the Fermi functions approximate delta functions, the above formula turns into : 



Ql(V) = -^—±[(A R + A c )ImG dd (eV) - XA R ImG M (XeV)] (3) 

7TftA 



and similarly: 



= d ^ = ~^ 1 ^lHA L +A c )ImG dd (XeV)-A L ImG dd (eV)} (4) 
av 7mA 

din 4e 2 An 

g c(V) = -i# = r-£[-A L ImG dd (eV) - XA R ImG dd (XeV)] (5) 

av 7mA 

The conductance Gl,r,c, defined in this way is a quantity that can be measured experimentally and helps for deter- 
mining the coupling ratios as we shall see below. 

We suggest now a way to search for Kondo dips in the differential conductance. One has to compute the second 
derivative of Ql,r,c for V = which contains the convexity information (in case the first derivative vanishes, the 
convexity determines if one has a peak or a dip). The second derivative of the differential conductance can be expressed 
as: 

Gl"(0) = ^^/mGL(0)[(l + ^) - A 3 ] (6) 
7mA i\ R 

Qn"(0) - -^^^M [(1 + ^)X 3 - 1] (7) 

Gc " {Q) = 4e '«hA AR ImG '^ 0) + A '] ; (8) 

In the above formulas — ImG dd (0) is the second derivative of the spectral function at equilibrium (i.e. for V = 0) and 
is always negative for T < Tk indicating the presence of the Kondo peak in the density of states. This property will 
be used to determine the peak or dip behavior of the conductance. 

The distribution (and sign) of the conductances in the three-lead device is controlled by the coupling to the ground 
lead Ac- In the mentioned experiments @, 0] , when the aim was to read the non-equilibrium DoS, Ac is very small 
compared to the other couplings; in this case the currents L and R have always opposite signs, the same being true for 
the second derivatives (6-7); the conclusion is that both conductances are maximum in module at zero bias. However, 
for large coupling Ac the above formulae indicate the possibility of a more complex behavior of the conductance. 

We identify first the case A = — 1 , as being the one commonly used, with the gate potentials symmetrically applied. 
In Eq.6, the straight parenthesis on the r.h.s is always positive for negative A, so one may conclude that Ql (0) has the 
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same sign as — 7mGj d (0), which is negative, indicating a peak also in the differential conductance Gl- For negative A, 
it turns out from Eqs.3,4 that Q L > and Qr < 0, meaning that the left and right currents flow in opposite directions. 
Then the above argument holds for the module of Qr which is also maximum at zero bias. 

Interestingly, for — An and A = —I, Qc is not zero as a consequence of the spectral asymmetry ImGdd{eV) ^ 
ImGdd(-eV) (see Eq.5, and also notice the different heights of the side peaks in Fig.l). The same non-zero value 
of the current through the grounded lead was discussed in (T3 |. for the sequential tunneling case, and attributed to 
the electron-hole asymmetry. Under the above mentioned conditions Qc (V — 0) = indicating a saddle point. In 
this case, the conductance through the central lead does not show a peak shape, being instead antisymmetrical about 
V = . 




Energy 



FIG. 1: The dot density of states for the case of symmetrically applied bias (A = —1). The three Kondo peaks correspond to 
the chemical potentials in the three leads. Notice that the side peaks are different in height. The parameters axe:V — 0.01, 
e d = -0.12, Ar = A L = 0.01, Ac = 0.02 A C /A L — 2, Ar — Al — 0.01. The temperature is T = TV /10, and T K = 3.2 • 10" 4 . 
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FIG. 2: Sketch of a quantum dot connected to three leads, with the bias asymmetrically applied on the left and right leads 
(V and AV, with A = 1/1.4/1.58/1.8/2.5). (a) The differential conductance through the left lead changes from peak to dip 
(around zero bias). The value of the asymmetry parameter A where the crossover happens is Ao = ^/l + Ac / Al = 1.58. (b,c) 
Differential conductance through the right and central lead show a Kondo peak behavior being maximum in module at V — 
The other parameters : Ar — Al = 0.075,A C /A L = 3, e d = -0.12, T = T K /W, and T K = 3.2 • 10~ 4 . 



A different situation occurs in the case A > 0. It is obvious that, without loss of generality, we may consider actually 
A > 1 meaning that, by convention, the smaller bias is applied on the left lead. One may now conclude from Eq.6 
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that the peak/dip behavior in the left lead is completely determined by the sign of (1 + Ac/A#) — A 3 . This result 
may be of potential interest for experimentalists. In particular, if one finds the crossover value Ao = yl + Ac/Ajj 
(when the peak changes into dip), the ratio Ac/A# can be immediately extracted. The mentioned coupling ratio is 
otherwise inaccessible to direct measurement. Such a peak-dip crossover can be seen in the Fig. 2a. The numerical 
plot (see below calculation details) indicates that the dip is considerably narrower than the Kondo temperature, and 
is situated in the middle of a larger peak. This means that even if the conductance through the left lead initially 
increases, it will eventually start to decrease for a larger bias. It is an expected behavior, since for large biases all 
Kondo correlations arc killed and all the conductances approach the (much lower) non-Kondo limit. 

The conductances through the right and central leads have the usual peak behavior (see Fig. 2b, c). One has to pay 
attention to the fact that the conductance through the central lead is negative due to the sign convention regarding 
the currents direction, but it is maximum in module at zero bias, and should be considered as a Kondo peak. 

We stress that the formulas Eqs.6-8 are quite general, and also hold outside the Kondo regime (for high temperatures 
or weak couplings), but in this case they are not very insightful. The variation of the differential conductance is very 
smooth outside the Kondo regime, and it is difficult to determine its convexity unambiguously. 



III. APPLICATION: NUMERICAL RESULTS FOR THE MULTI-LEAD ANDERSON MODEL 

In order to illustrate the results, we shall calculate the differential conductance for the multi-lead Anderson model: 
H = Y^ zd&\&c + Ud\dfd\di + ^2 tk*cl aa c ka<T + ^ V a {c\ a(y d a + H.c), (9) 

where ea is the energy of the dot level, U is the Coulomb repulsion and, for instance, c kaiJ annihilates an electron 
with momentum k and spin a in the a lead. The last term describes the coupling between the dot and leads. The 
quantity needed in the current formula Eq.l is the dot retarded Green function Gdd,a =<< d^d], >>, which will 
be calculated in the Lacroix 0] approach. We write the equation of motion for the Green function and perform a 
decoupling approximation at the third step in order to close the system of equations. In our case, wc consider three 
leads (reservoirs) instead of one as in 0], the generalization being straightforward. The resulting formula, in the 
U — y oo limit, reads: 

G dd , a {uj + iQ) = ? ^ (10) 

, , w _ v 1 -AH — v v v < v-° Ct ' i -° > _i_ _AS_ v v < d -* c <"*-°> 

w a l^k.a u)-t k 2-ik,q,a,P V <* V P w-e k ^ ^k,a uj-e k 2^k,a V <* uj-e k 

Since we work in a non-magnetic case, the subscript a can be omitted in the Green function notation. It comes now 
into question how to treat the averages of the type < <y q(x _ a Ckp- a >, for instance, where the two operators refer to 
electrons in different leads. For equilibrium, when all the Fermi levels are equal, one can use the fluctuation-dissipation 
theorem which relates the average to the corresponding Green function: 

< c\ a -o c kp-<7 >= ~- J f(u)Im « c kfj - a c\ a _ a » duj, (11) 

In non-equilibrium, the Fermi functions f a arc different in the three leads a = L, R, C, and for the calculation of the 
mentioned average we propose 

< c\ a _ a c k ^ a >= -(1/2tt) J (f a (w) + fp(u))Im « c kl 3^c\ a _ a » duj. (12) 

In the calculations there also appear averages of the type < d}_ a c ka ^ a >, so that a distribution function on the dot 
is also needed. For the dot we shall use the weighted mean fd = J2 a A Q / Q /A. This intuitive expression of /<j has 
been used in, e.g. (l3l. [l4|. and we have checked that the present results, i.e. existence of the dip and its temperature 
dependance, are stable with respect to different ways to introduce the non-equilibrium. 

These non-equilibrium aspects being settled, we continue along the Lacroix line, and with the definitions A = 
£A Q = jV = „K 2 , A a = f G y )J °("'W , and B a = ^ f i^\ du' (where D is the 

half-bandwidth and is taken as energy unit) the following expression is obtained : 

r ( l-^ + GS> + iO)Eafr» ri „. 
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In getting this expression we use the approximation A a w ~G* dd B a , which was proposed first by Lacroix for the 
Kondo regime, and later used also in 

[U Eg. E q 13 

can be easily solved. First, one conjugates the equation, and after 
that G* d (to + iO) is replaced in the initial formula to obtain a simple algebraic second degree equation for Gd(uJ + iO), 
only one of the solutions being physical. The solution is analytic, but is not included here because its lengthy. It is 
important to mention that this solution of the EOM allows to address the low temperature regime (T <C Tk), which 
is not possible with the more simple decoupling scheme in [3| ■ 
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FIG. 3: Evolution of the Kondo dip with temperature; A = 2; the other parameters as in Fig. 2 

Fig. 3 shows that the Kondo-type temperature dependance of the conductance (i.e., decrease of with increasing 
temperature) occurs in the range eV/Tjc € [—0.5, 0.5] (Tk being calculated as in This range is in fact dependent 

on A and we checked that it becomes narrower when A increases. 




FIG. 4: The dip width scaled with temperature, as function of the asymmetry parameter A, for three different temperatures 
below Tk- The linear temperature dependence is obvious for A < 2.4, where the three scaled curves coincide. The dip vanishes 
at Aq = 1.58 which can be obtained by extrapolation. 



One can also notice that the width of the dip increases with the temperature. Actually the temperature dependance 
of the conductance comes from the spectral function ImGdd(V = 0,T). Let us define the width of the dip w(X,T) 
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as the distance between the two peaks. By scaling the width with temperature the data from Fig. 3 generates Fig. 4 
which shows that, for not too large values of the parameter A, the width behaves like: 



w(X, T) = /(A) T, with /(A = A ) = 0. (14) 

For A close to Ao the error in the calculation of the dip width becomes large, and we suggest a numerical fit to obtain 
the dip-peak crossover point Aq. The polynomial fit gives a value of Ao very close to the exact result Ao = 1.58. 



IV. CONCLUSIONS 



In conclusion, we have studied the non-equilibrium transport through a quantum dot connected to three reservoirs. 
In Kondo regime, if the potentials are applied asymmetrically on the leads (i.e. Vr = AVj,, Vc = 0, A > 1) a peak-dip 
crossover of the differential conductance occurs in left lead at zero bias. The proof is based on the connection between 
the convexity of the differential conductance and the convexity of the density of states at V — and is independent 
of the Hamiltonian model. The ratios of the lead-dot couplings (A^/Ac, Ajj/Ap) can be obtained if the crossover 
value Aq is determined experimentally. 
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